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1 .  Introduction 


This  research  considers  only  passive  vibration  isolation.  Most 
physical  structures  designed  for  dynamic  environments  have  isolator 
elements  to  attenuate  the  response.  Examples  of  problems  that  could 
benefit  from  this  research  are  ground  vehicle  response  and  equipment 
or  instrument  vibrational  response.  These  problems  use  passive  isola¬ 
tors  to  minimize  the  vibrational  response. 

As  a  preliminary  step  toward  the  study  of  active  structural  vibra¬ 
tion,  the  passive  vibration  is  studied  in  this  research.  Most  Aerospace 
structures  of  the  future  will  use  active  vibration  isolation  to  attenuate 
the  response.  The  active  vibration  isolation  problem  would  contain  the 
same  constraints  used  in  this  study  with  the  addition  of  control  forces 
to  the  state  equations  and  to  the  objective  function  -  Large  space 
structures  would  use  active  vibration  isolation  to  control  the  response 
of  the  very  flexible  structural  systems. 

The  algorithms  that  are  developed  in  this  research  are  consistent 
with  the  minimum  weight  problem  which  is  not  considered  in  the  present 
study. 

This  research  considers  constraints  that  are  displacements, 
accelerations  and  natural  frequencies.  The  design  variables  are  linear 
changes  to  mass,  stiffness  or  damping  matrices.  The  constraints  can 
be  expressed  in  either  the  time  or  frequency  domain  and  the  cumulative 
constraint  is  used  to  measure  the  amount  of  constraint  violation.  It 
is  shown  that  the  variation  of  the  displacements  or  acceleration  constraints 
are  shallow  in  reciprocal  design  variables.  The  objective  function 
represents  a  design  variable  that  restrains  displacements  or  accelerations 
to  be  less  than  a  maximum  value. 

These  algorithms  have  been  studied  for  transient  response,  frequency 
response  and  stationary  random.  No  attempt  was  made  to  consider  ill 
conditioned  vibrational  problems  that  occur  between  a  structure  and 
isolator  elements  that  are  designed  to  move  independent  of  the  structure. 


2.  Transient "Response 

The  minimization  of  displacements  or  accelerations  can  be 
formulated  as  a  MIN-MAX  optimization  problem. 
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Only  the  direct  method  of  solution  has  been  considered  in  this 
study.  These  algorithms  can  be  used  with  the  modal  formulation  and 
a  separate  section  references  the  recent  work  completed  using  modal 
analysis. 

Equation  (1)  minimizes  the  maximum  acceleration  in  the  time 
domain.  The  objective  function  could  be  displacements  instead  and  the 
present  algorithms  could  be  used  directly.  Equation  (2)  is  the 

eigenvalue  problem  whose  inclusion  permits  the  use  of  frequencies  in 


equation  (5)  in  the  analysis.  Equation  (3)  is  the  structural  dynamic 
equations  in  matrix  form  which  describe  the  displacement  response  X(t) . 
Equation  (4)  is  the  so  called  relative  displacement  or  rattlespace 
constraint.  The  present  algorithms  can  include  this  type  of  constraint 
in  the  analysis.  However,  no  specific  numerical  examples  are  presented 
using  the  rattlespace  constraint.  Equations  (6),  (7),  and  (8)  show 
the  linear  changes  to  the  stiffness ,  mass  or  viscous  damping  matrix 
with  the  design  variables  cu.  The  design  variables  could  contain 
differing  sets  in  equations  (6) ,  (7) ,  (8) .  Equation  (9)  lists  the 
constraint  limits  on  the  design  variables  a. . 


Sometimes,  it  is  convenient  to  solve  vibration  problems  in  the 
driving  frequency  (a  domain.  This  is  true  for  problems  which  have 
experimentally  available  results  for  transfer  functions.  Also,  for 
stationary  random  analysis,  the  frequency  domain  transfer  function 
must  be  determined.  Equation  (3)  is  transformed  to  the  steady 
state  frequency  domain  by, 

(10)  X  «  RE  {X  e^1}  ,  P  -  RE  {  P  e^1} 

o  o 

where  KE:  denotes  real  part  of 
i  ■  /-I 

<•>:  driving  frequency 
X  :  amplitude  of  harmonic  response 
Pq:  amplitude  of  harmonic  loading 

For  the  harmonic  substitution,  equation  (3)  becomes, 

(11)  (-u2**  +  it»)C  ♦  K)Xq  -  po 

The  amplitudes  Xq,  Pq  axe  complex  numbers.  Equation  (11)  may  be  solved 

repeatedly  for  Xq  given  Pq  and  oi  using  complex  arithmetic.  It  is  more 

convenient  to  use  the  real  displacement  components  in  the  analysis.  The  method  of 

reference  (1)  is  used  to  work  with  the  real  and  imaginary  components 

of  X  . 
o 


Hie  optimization  becomes  in  the  frequency  domain. 
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Equation  (13)  is  the  amplitude  of  steady  state  acceleration  and  equation 
(15)  axe  the  structural  dynamic  equations  to  be  solved. 


4.  Stationary  -Random 

A  Frequency  Response  solution  is  first  analyzed  to  determine  the 
transfer  function  H(co)  which  is  either  the  displacement  or  acceleration 
at  a  response  point  of  interest.  The  spectral  density  of  the  output 
is  given  in  terms  of  the  spectral  density  of  the  input  for  a  single 
input/output  system  is  given  in  reference  (2) • 

SQ(w)  =  |  H(a>)  |  2SjG«0 

The  same  reference  also  lists  techniques  for  analyzing  multiple 
input/output  systems.  The  mean  square  value  can  be  calculated  for 
any  frequency  interval, 

W2 

-2, 
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U1 

Various  performance  measures  have  been  proposed  for  random  analysis 
such  as  using  either  the  spectral  density  or  mean  square  value.  The 
optimization  problem  for  stationary  random  becomes. 
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The  maximum  displacement  or  acceleration  spectral  density  is  the 
objective  function  to  be  minimized  in  equation  (22).  Only  the  single 
input/output  case  is  used  to  calculate  the  spectral  density  by  equation 
(25),  the  objective  function  is  converted  to  a  set  of  equivalent 
integral  constraints  and  the  minimization  is  done  then  on  the  mean 
square  response  in  effect. 


5.  Cumulative  Constraint  and  MIN-MAX  Problem 


The  cumulative  or  equivalent  integral  constraint  has  been  used  in 
the  optimal  control  literature  (  3  )  to  convert  many  discrete  points  in 
the  time  domain  to  one  equivalent  integral.  Thus  many  discrete  constraint 
equations  are  lumped  into  one  equation.  The  cumulative  constraint  measures 
the  total  amount  of  constraint  violation.  If  a  satisfied  constraint  is 
of  the  form, 

\ 

(31)  <J>(t)  <0  all  t. 

The  brae  function  measures  the  amount  of  constraint  violation, 

<|>(t),  4>(t)>  0 

«Kt)>  »  { 

0  ,  $(t)<  0 


A  constraint  totally  equivalent  to  (31)  is, 

(32)  /«Kt)>  dt  *  0 

Equation  (32)  is  not  identically  zero  if  the  constraint  (31)  is 
violated.  Instead  of  using  (31)  at  many  discrete  points  in 
time,  one  total  constraint  (32)  is  used  which  measures  where  the 
constraint  (31)  is  violated. 

The  objective  function  (1),  (13)  or  (22)  can  be  converted  to  a 
simplier  algebraic  form. 

Consider  equation  (1), 


MIN  (MAX  |X(t)|  ) 

This  minimization  is  equivalent  to  minimizing  an  additional  design 
variable  a  such  that 
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MIN  a 


|x(t)|—  a<  0  for  all  t 

The  acceleration  constraint  that  was  just  introduced  can  be  converted 
to  a  cumulative  constraint. 

The  MIN-MAX  part  of  the  optimization  becomes. 


MIN  a 


/<  |X(t]|  -  a  >  dt  -  0 


To  account  for  the  absolute  value  of  acceleration,  the  integral 


is  written 


fl  dt  *  0 


I«X-a,X>a 

•*  •• 

I  =  -X-  a ,  -  X  >  a 

I  *  0  otherwise 

The  objective  functions  (13)  or  (22)  in  the  frequency  domain  are 
computed  in  the  same  manner  with  frequency  replacing  time  in  the  integral. 

The  inner  problem  or  the  maximization  in  this  research  was  done 
by  function  evaluation.  This  is  efficient  for  the  transient  problem, 
but  the  frequency  response  problem  requires  a  decomposition  for  each 
driving  frequency  in  equation  (IS) .  It  would  be  required  to  reduce  the 
basis  of  equation  (15)  by  using  the  real  normal  modes  for  efficient 
solution  in  locating  the  maximum.  Reference  (1)  recommends  performing 
a  one  dimensional  search  on  the  variable  u  to  locate  the  maximum.  This 


S 


6.  Reciprocal  Design  Variables 


For  statically  determinate  structures ,  stresses  and  deflections 
are  proportional  to  design  variables  that  are  linear  changes  to 
stiffness  such  as  areas  of  rods  in  truss  members.  For  indeterminate 
structures,  this  is  only  an  approximation.  It  was  investigated  in 
references  (4,5)  and  found  that  high  quality  explicit  expressions 
for  stresses  and  deflections  could  be  generated  using  a  first  order 
Taylor  series  expansion  in  reciprocal  design  variables.  That  is, 
the  design  variable  space  for  stress  and  static  deflection  is  shallow 
in  reciprocal  design  variables.  The  linearized  Taylor  series  expansions 
represent  lines  that  are  very  good  approximations  to  the  exact  constraints. 
The  expansion  of  a  response  quantity  $  is  done  in  the  reciprocal  design 
variables  &  , 
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The  mass,  stiffness  and  damping  matrices  are  linear  functions  of 
the  design  variables  a. .  The  derivative  of  a  general  matrix  G  in  the 
above  classification  is, 
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The  derivative  can  be  calculated  by  setting  the  linear  property 
equal  to  one  in  the  element  matrix  G,  and  multiply  by  the  square  of 
the  design  variable.  The  same  subroutine  used  to  calculate  element 
matrices  can  be  used  to  form  the  above  derivative. 


.  Analytic  Derivatives 


The  direct  solution  of  the  dynamic  response  equations  in  the 
time  domain  uses  an  efficient  implicit  equation  solver  such  as  Newmark 
integration.  The  same  technique  could  be  used  in  the  modal  formulation 
if  coupled  real  mode  damping  matrix  exists.  This  would  be  the  most 
general  capability  for  solution  of  the  dynamic  equations.  The  Newmark 
integration  equations  (6)  are  listed  for  one  set  on  integration 
parameters. 


h,  h, 

(4G) 

MX  +  CX  ♦  KX  p 

Z2  Z2  Z2  ‘  *t2 

(41) 

KK  *  K  M 

+  (feC 

(42) 

PP„  »  P^  +mL4ts2  X^  +  X.  +  X  1 

t2  t2  |fAt)  tj  St  tj  tjj 

+  C 

ks  *  \  1 

(43) 

KK  *X.  -  PP. 

Z2  2 

(44) 

\  mm2  c  x,2  ■ 

(45) 

Y  Y  At  "  At  L' 

xt,  xt,  "J  xt,  *  1  xt. 

The  displacements  at  the  next  time  step  are  calculated  by  equation 
(43) .  The  matrix  KK  is  only  factored  when  the  time  step  At  changes . 

The  acceleration  and  displacements  are  recovered  by  equations  (44)  and 
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(45).  The  derivatives  of  the  response  quantities  (7)  are  found  by 
differentiating  equations  (47)  THRU  (45) .  This  is  the  pseudo  loads 
technique.  The  derivatives  with  respect  to  the  reciprocal  variables 
are, 

fA ~  ax.  app* 

ae.  *bA2 

(47)  3KK  3K  4  ,  3H  2  3C 
367  *  36.+  m2  W:  *6. 


C49)  3X 

—2  *  4  , 

36*  TS~t)2 


(50)  3X  3X  jfX. 

— 2.— 1  +  A t  — 1 
36i  36*  2  36i 


•• 


Using  this  technique,  the  derivatives  of  displacement,  velocity 
and  acceleration  must  be  calculated  and  saved  for  all  degrees  of  freedon 
in  the  finite  element  model  at  two  neighboring  points  in  time.  The  KK 
matrix  in  (46)  was  decomposed  in  the  response  calculations  and  would  not 
be  factored  again  in  this  step. 

The  pseudo  loads  technique  was  applied  to  the  structural  equations 
(24)  in  the  frequency  domain. 
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As  discussed  in  reference  (8)  ,  frequency  constraints  are  inherently 
nonlinear,  Reference  (9)  suggests  if  the  structure  has  large  fixed  masses, 
direct  design  variables  and  a  linear  Taylor  series  expansion  should  be 
used.  To  use  frequency  constraints  in  the  present  research,  the  inverse 
eigenvalue  A*  should  be  used, 

KY  =  AMY 

A*KY  »  MY 


The  justification  follows  from  the  Rayleigh  quotient: 


YTM  Y  +  Za.YTM.Y 
A*  =.  - i—i_ 

YXK  Y  +  &*.YTK.Y 
o  11 


where  the  mass  and  stiffness  are  given  by  equations  (6) ,  (7)  or 
(18),  (19). 

For  large  fixed  masses,  M  dominates  the  top  of  the  ratio.  The 
reciprocal  eigenvalue  contains  cl  on  the  bottom.  The  eigenvectors  Y 
are  functions  of  the  design  variables,  but  enter  in  the  same  manner 
at  the  numerator  and  denominator.  A*  might  be  best  represented  in  0. 
space. 

The  derivative  is  well  known  and  given  in  reference  (10) : 


(52)  aAi 

3oj 


pK 

ji«r 


-  uiSL 


This  derivative  has  to  be  expressed  in  terms  of  the  reciprocal 
design  variables  and  reciprocal  eigenvalue. 
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Linearized  Constraints 


The  acceleration  cumulative  constraint  has  the  following  first 
order  Taylor  series  expansion  in  the  reciprocal  variables. 


/T  I  dt  =  0 
o 


where  I  =  X  *  T-It—  6$.  -  6»  for  X  -  g 

O  1 

I  *  -X0  --  I-§~-  36  -  6,  for  -X  >  2 

1=0  For  X  otherwise 


This  constraint  is  numerically  integrated  by  a  modified  trapezoid 
law  which  finds  those  response  points  above  a  line  for  which  the  constraints 
are  violated.  The  algorithm  interpolates  to  find  the  points  where  the 
actual  constraint  is  violated. 

The  steady  state  acceleration  amplitude  in  the  frequency  domain  is, 


A  =  oj2  /U2  ♦  V2 


where 


XQ  » U  -  i y  is  the  amplitude  of  steady  state  displacement . 


The  first  order  Taylor  series  expansion  for  the  acceleration 
amplitude  magnitude  is. 


where 
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The  acceleration  spectral  density  of  the  output  can  be  expressed 
in  terms  of  the  spectral  density  of  the  input  and  the  transfer  function 
for  a  single  input/output  system. 

SQ  (to)  ■  a)4  |  H(w)  j  2Sj(u) 

Where  the  transfer  function  is  the  displacement  at  the  response  point. 

To  utilize  the  same  approximation  as  the  frequency  response 
acceleration  constraint  equation  (55) ,  consider  the  square  root  spectral 
density. 


S«  JSq  (to)  *w2  |H  (o»)| 

2 

To  modify  (55)  for  acceleration  spectral  densities,  the  to  is 
multiplied  by  JS'j  (to) . 


(56)  j  du  =  o 


where 


I  *  S_  + 


to  ✓  Sj  (to) 


✓  IT  ♦  V 


A  formula  could  be  derived  for  multiple  sources  with  cross  correla¬ 
tion  in  a  similar  manner.  Equation  (56)  developed  with  this  algorithm 
is  the  square  root  mean  square  acceleration  that  exceeds  a  specified 
value  in  the  frequency  domain. 
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9.  Modal  Method 


The  direct  method  was  analyzed  in  this  research,  but  the  developed 
algorithms  could  be  used  with  the  modal  analysis  method  with  the  addition 
of  the  following  algorithms. 

The  structural  dynamic  equations  are  first  converted  to  modal 
amplitudes  as  the  degrees  of  freedom. 

X  «  YZ 

•• 

(57)  (YTMY)Z  +  (YTCY)Z  +  (YTKY)Z  «  YTP 

The  generalized  mass  and  stiffness  are  diagonal,  but  the  modal 
damping  is  in  general  a  coupled  matrix  as  the  result  of  being  transform 
by  the  real  eigenvectors  Y.  When  derivatives  are  calculated  for  equation 
(57)  it  becomes  necessary  to  differentiate  the  eigenvector.  The  research 
presented  in  reference  (11)  should  be  used  to  efficiently  calculate  the 
eigenvector  derivative  if  large  changes  have  been  made  to  the  structure. 
If  only  isolator  changes  are  made,  usually  only  small  changes  are  made 
to  the  structure. 

In  such  a  case,  the  eigenvector  derivative  would  not  be  needed. 

A  direct  modal  solution  was  done  for  a  problem  of  this  type  in  reference 
(12).  In  this  problem^ the  information  could  be  accurately  calculated 
in  the  Ritz  eigenspace  of  the  current  design.  A  direct  modal  solution 
could  also  be  used  with  the  present  research.  The  algorithms  of  the 
present  research  should  replace  the  very  inefficient  one  developed  in 
reference  ( 12) . 


10.  Sequential  Linear  Programming 

The  problem  considered  in  this  research  of  minimizing  a  linear 
design  variable  subject  to  constraints  on  displacements  or  accelerations 
in  the  time  or  frequency  domain  is  an  almost  linear  problem  in  reciprocal 
design  space.  It  is  only  natural  to  use  sequential  linear  programming 
as  the  optimization  algorithm.  A  primal-dual  linear  program  which  is 
listed  in  reference  (13)  was  used  as  the  optimizer.  Sequential  linear 
programming  is  described  in  reference  (14).  Reference  (15)  has  shown 
the  dynamic  response  optimization  problem  to  have  a  disjoint  design 
space.  Sequential  linear  programming  (14)  is  capable  of  solving  such 
a  nonconvex  problem. 

The  design  variables  (9),  (21)  are  converted  to  the  reciprocal  space. 


(58) 


The  actual  design  variables  were  the  changes  $0^, 
(59)  8i  -  \  ^ 

a. 


where  "0.  is  the  current  reciprocal  design  variable. 

The  design  variables  90^  are  not  restricted  in  sign.  Linear 
programming  requires  the  design  variables  to  be  non-negative.  This 
requires  the  design  variable  to  be  converted  into  the  difference  of 
two  non-negative  variables. 

(60)  5^  • 

>  0,  >0 

The  bounds  on  $5^  become, 
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(61) 


1 


aiu  a. 

1 


<  66i  -  5^  <  1 _  .  1_ 

aiL  a. 


To  maintain  linear  approximations,  the  design  variable  is  also 
bounded  by  a  move  limit  NL  which  is  taken  as  a  percent  change  from  the 
current  design  point. 


60.  <  M. 

i  “  i 


The  bounds  on  the  new  design  variables  in  equation  (61)  become. 


0  <  60t  <  MIN  M  ,  (-^ - hr  ) 


iL  a. 
i 


0  <  667  <  MIN  M. ,  (-  i—  +  -i— ) 


iu  a. 
x 


The  remainder  of  the  optimization  algorithm  is  converted  by 
equation  (60)  to  the  66?,  66 7  space.  For  example,  the  acceleration 
constraint  (54)  in  the  time  domain  becomes. 


where 


/  I  dt  *  0 


1  *  Xo  *  2  (  66i  • 

i 


H6i) 


for  x  >  6 


1  <  5Bi  -  H?  e 


for  X  otherwise 


-X  >  6 


11.  Ill  Conditioning 


As  mentioned  previously,  the  vibration  isolation  problem  can  be 
ill  conditioned.  Isolators  are  designed  to  move  independently  of  the 
main  structure.  This  would  require  the  properties  to  be  several  orders 
of  magnitude  less  than  the  surrounding  structure  and  conditioning 
problems  would  then  arise.  An  approximate  way  of  dealing  with  this 
situation  which  was  done  in  reference  (12)  is  to  use  modal  synthesis . 
The  main  structure  is  modeled  by  component  modes  and  the  isolators 
are  included  at  the  syntheses  time.  If  the  generalized  mass,  stiffness 
and  damping  are  of  the  same  order  of  magnitude  as  the  isolator  proper¬ 
ties,  the  problem  will  become  well  conditioned.  However,  this  method 
is  approximate  and  no  error  estimates  are  available  for  it. 

Large  space  structures  are  very  flexible  and  inherently  ill 
conditioned.  More  accurate  techniques  would  have  to  be  used  for  prob¬ 
lems  in  this  category, 

A  method  for  solving  ill  conditioned  eigensysteans  is  presented 
in  reference  (16) ,  but  it  may  require  modifications  to  be  efficient 
for  structural  problems.  A  method  for  the  solution  of  ill  conditioned 
stiffness  matrices  is  presented  in  reference  (17)  and  it  might  be 
possible  to  modify  it  for  the  generalized  eigenvalue  problem.  The 
numerical  integration  schemes  for  stiff  equations  are  presented  in 
reference  (18) . 


The  model  of  reference  (1)  shown  in  figure  1  was  subjected  to  the 
displacement  inputs  f^(t)  and  fg(t)  shown  in  figure  2.  This  model 
represents  a  vehicle  running  over  a  bump.  The  transient  step  size 
used  was  .1  sec  and  39  time  intervals  were  calculated.  The  five 
springs  were  used  as  design  variables  with  limits  shown  on  figure  1. 
The  objective  function  was  a  design  variable  which  represented  the 
maximum  acceleration  at  point  1  in  the  model  over  the  3.9  sec  time 
of  response.  The  acceleration  constraints  were  made  active  when 
the  acceleration  bound  was  99%  of  the  maximum.  Figure  3  presents 
the  decrease  in  acceleration  at  point  1  in  the  model  versus  the 
required  number  of  structural  analyses.  By  equation  (60) , each  design 
variable  is  converted  into  two  so  the  total  number  of  design  variables 
used  in  the  linear  program  was  eleven.  Initially,  the  reciprocal 
variables  were  constrained  by  a  move  limit  to  lie  within  *25%  of 
the  initial  values.  Convergence  was  obtained  at  iteration  three. 

The  spring  rates  found  at  the  optimum  were, 


k2  »  51.2  lb/ in 
kg  *  200,  lb/in 
kj  -  200  .lb/in 
k4  «  1600.  lb/in 
k5  »  1000.  lb/in 

2 

The  minimum  acceleration  obtained  was  228.8  in/sec  .  Figure  4 
presents  the  initial  response  versus  the  optimal  one. 
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An  indeterminate  beam  composed,  of  rods,  shear  panels,  springs 
and  concentrated  masses  is  shown  in  figure  5.  The  two  point  masses  were 
subjected  to  a  harmonic  force  input  of 

1000  SIN  (  ir  t)  0  <t  <  sec 

f(t)  =  { 

0  2  sec  <t 

Eight  design  variables  were  used  in  the  analysis,  3  rod  areas 
at  top  and  bottom  and  the  two  scalar  springs.  The  two  vertical  rods 
inside  the  shear  wall  remained  unchanged  at  the  initial  area  of  the 
assembly.  The  vertical  displacement  response  of  the  top  two  response 
points  were  required  to  be  less  than,  the  9th  design  variable.  This 
problem  validated  the  extension  of  the  method  to  finite  elements  and 
multiple  response  points. 

This  problem  is  also  statically  indeterminate  and  the  reciprocal 
design  variables  used  in  the  analysis  were  only  an  approximation,  A 
lumped  mass  model  was  used  for  the  rods  and  shear  panels.  The  shear 
panel  used  was  the  hybird  one  derived  in  reference  (19),  The  damping 
used  was  C  *  .01  LB-SEC/IN  and  a  scalar  viscous  damper  of  this  magnitude 
to  ground  was  applied  at  every  degree  of  freedom.  Figure  6  presents 
the  decrease  in  the  maximum  displacement  vs  structural  iteration. 
Convergence  was  achieved  at  iteration  3  when  the  maximum  displacement 
went  from  an  initial  value  of  .0166  IN  to  .0082  IN  at  the  third  iteration. 
The  rods  at  the  optimum  all  reached  the  upper  bound  of  2  IN  and  the 
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scalar  springs  had  the  value  of  51,2(10  LB/IN  which  were  near  the  lower 
bound  but  not  exactly  equal  to  it. 

Throughout  the  optimization,  the  symmetry  in  the  problem  was 
preserved.  This  problem  could  be  ill  conditioned  when  the  spring  rates 
were  more  than  a  power  of  ten  smaller  than  the  lower  limit.  Typically 
isolator  elements  are  designed  to  vibrate  independently  from  the  surround¬ 
ing  structure.  Vibration  isolation  is  inherently  ill  conditioned  for  this 
reason*  and  this  research  has  not  considered  the  ill  conditioned  problem. 


Frequency  Response: 


The  model  shown  in  figure  1  was  subjected  to  equal  inphase  displace^ 
ment  inputs  fj(t)  *  f2(t)  *  5  coswt  at  the  tire.  This  would  represent  a 
vehicle  on  a  shaker  table.  The  five  springs  were  used  as  design  variables 
with  the  limits  shown  of  figure  1.  The  acceleration  amplitudes  were 
evaluated  between  5  RAD/SEC  and  44  RAD/SEC  in  steps  of  1  RAD/SEC  as  the 
driving  frequencies.  The  objective  function  was  the  maximum  steady  state 
acceleration  amplitude  at  point  1  over  the  range  of  driving  frequencies. 

The  acceleration  constraints  were  made  active  when  the  acceleration  was 
99%  of  the  maximum.  Figure  7  presents  the  decrease  in  acceleration 
amplitude  at  point  1  in  the  model  versus  the  required  number  of  structural 
analyses.  Initially,  the  reciprocal  variables  were  constrained  by  a  move 
limit  to  lie  within! 25%  of  the  initial  values.  When  a  step  was  not 
minimizing,  the  percent  move  limit  on  the  reciprocal  variables  was 
decreased  by  50%  and  the  linear  program  was  resolved  at  the  previous  design 
point.  Convergence  was  achieved  at  iteration  7  which  was  close  to  the 
value  found  at  iteration  5. 

The  spring  rates  at  the  optimum  were  as  follows: 

kj  •  52.9  LB/IN 
k2  *  231.2  LB/IN 
k3  -  215.1  LB/IN 
k4  -  1000,0  LB/IN 
*5  *  1311,0  LB/IN 
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Stationary  Random  Response: 


The  model  shown  in  figure  1  was  subjected  to  a  random  displacement 
at  the  tire  patches  as  discussed  in  reference  (20)  with  parameters  that 
correspond  to  a  smooth  highway.  A  Frequency  Response  solution  is  first 
completed  with  a  unit  harmonic  displacement  elb*  at  the  front  tire  and 
a  unit  displacement  with  phase  lag  e*  at  the  rear  tire  with 

phase  angle  <f>  ~  .  The  spectral  density  of  the  output  in  terms  of  the 

spectral  density  of  the  input  and  transfer  function  is, 

So(a»)  -  |H(o>)  I2  SjCotf. 

The  transfer  function  is  determined  by  using  the  acceleration  output 
of  the  frequency  response  solution  due  to  the  unit  harmonic  input. 

The  spectral  density  acceleration  was  evaluated  between  5  RAD/SEC  and 
44  RAD/SEC  in  steps  of  1  RAO/SEC.  The  objective  function  was  the  design 
variable  representing  the  maximum  acceleration  spectral  density. 

The  acceleration  constraint  was  made  active  when  it  was  99%  of  the 
maximum.  Figure  9  presents  the  decrease  in  the  objective  function  versus 
the  required  number  of  structural  analyses.  The  reciprocal  spring  rates 
were  constrained  by  *  25%  of  the  current  value  as  move  limits.  At  the 
detection  of  each  infeasibility,  the  move  limit  was  reduced  by  50%  of 
the  current  percent. 

Convergence  was  achieved  at  iteration  3. 

The  spring  rates  at  the  optimum  were, 


Kj  =  51.2  LB/IN 

K2  »  200.  LB/IN 

K3  =  200.  LB/ IN 

K.  »  1000.  LB/IN 
4 

Kg  =  2000.  LB/IN 

2  2 

The  minimum  spectral  density  was  49.61  (IN/SEC  )  /Hz.  The  initial 


and  optimized  spectral  densities  are  presented  in  figure  10. 


Transient  Response  with  Frequency  Constraints: 

The  model  shown  in  figure  1  was  again  used  with  the  same  constraints 
and  loading  condition.  An  additional  constraint  was  added  on  the  inverse 
frequency: 

KY  -  X*MY 


A  10%  move  limit  was  imposed  initially  on  the  design  variables  and 

each  time  an  infeasibility  was  found,  the  move  limit  was  decreased  by  50%. 

2 

The  minimum  acceleration  found  was  286.9  IN/SEC  after  six  structural 
analyses.  The  analysis  was  done  in  both  the  direct  eigenvalue  and  inverse 
eigenvalue  space  with  the  same  convergence  properties.  Figure  11 

presents  the  decrease  in  acceleration  vs  structural  analysis  in  the  inverse 
eigenvalue  constraint.  Near  convergence  was  achieved  at  four  analyses, 
but  no  conclusions  concerning  the  eigenvalue  constraint  could  be  inferred 

from  this  problem. 


13.  Conclusions 


The  vibration  isolation  problem  can  be  efficiently  analyzed 
expressing  the  displacement  or  acceleration  constraints  as  a  first  order 
Taylor  series  expansion  in  reciprocal  design  variables.  The  cummlative 
constraint  folds  many  discrete  constraints  into  a  single  one. 

The  method  for  approximating  the  displacement  or  acceleration 
constraints  is  consistent  with  the  current  algorithms  for  weight 
minimization. 
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M^g  =  290  ib. 

M^G  =  4,500  ib. 

I  =  41,000  lb. -in- sec 
tyG  =  96.6  lb. 

r^G  =  96.6  lb. 


Cx  =  10  lb. -sec/in 
-  25  lb.-sec/in 
Cj  s  25  lb.-sec/in 
Qj  *  5  lb.-sec/in 
Cij  —  5  lb.-sec/in 


LOWER  LIMITS 
W-2  =  50  lb. /in 

=  200  lb. /in 
=  200  lb. /in 
a  MX)  lb. /in 
=  1000  lb. /in 


INITIAL  DESIGN 

*  100  lb. /in 

l<2  *  300  lb. /in 
-  300  lb. /in 
=  1500  lb. /in 

*  1500  lb. /in 


UPPER  LIMITS 
KU^  =  500  lb. /in 

KU^  =  1000  lb. /in 

KUj  =  1000  lb. /in 

KUjj  =  2000  lb. /in 

KUij  =  2000  lb. /in 


FIGURE  1:  OPTIMIZATION  MODEL 
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AMPLITUDE  XQ  =  5* 

VEHICLE  SPEED  S  =  450  IN/SEC 

WHEEL  BASE  L  =  120  IN 

=  360",  d2  =  144" 
w,  =  jrs  -  1.25  T r 

w-  -  jffS  =  3.125TT 

*i  =  div  ,3  SEC'  ta*  (d’+d2>s  *1.12  SEC 

TIME  LAO  FRONT  TO  REAR  tL 
*L  =  u/$  =  .2667  S€C 

FRONT  WHEEL  DISPLACEMENT 

ft  (t)  =  X0(l  -  COSWt)  0*Utf 

!,(»)  =  X0(l  +  COSW(t-tj))  tt*t*t2 

REAR  WHEEL  DISPLACEMENT 

f2(t)  =  f,  (t-ta)  <Kt-tL*t2 

FIGURE  2:  TRANSIENT  DISPLACEMENT  INPUT  FOR  FIGURE  1 
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i 


LOWER  LIMIT 
.5  IN2 


INITIAL  AREAS 
1  IN2 


E  = 

10x10®  ps I 

M  = 

G  = 

4xl06  ps  I 

f  - 

8.79x  10'3  SLUGS/IN  3 

C  = 

t  = 

.1  IN 

LOWER  LIMIT  INITIAL  K 

50,000  LB/lN,  100,000  LB/IN. 

FIGURE  5:  RODS,  SHEAR  PANELS,  SPRINGS  AND  MASSES 


UPPER  LIMIT 
2  IN2 

.005  SLUGS 

.01  LB-SEC/IN 

UPPER  LIMIT 
200,000  LB/IN. 
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